Author: Ranjit Deshmukh 11/13/2016 This script does the following 1) Creates the results table using outputs from the results csv from the Renewable Energy Value model
## [1] "/Users/ranjitster/Dropbox/renewable_energy_value/renewable_energy_value"
Input files
Specify scenarios to process
User parameter inputs
Function to grab legend
Import results and other parameter from csv files
Creating empty results data tables and lists for plots.
Capital recovery factor
Preparing the results input table. Splitting the scenario and scenario_econ_dispatch names into 2. Replacing the No Re scenario builds S0W0 for battery storage scenarios with the base S0W0 scenario numbers.
VRE is variable renewable energy.
total_cost_new_conv: Total annual cost of new conventional capacity (annualized cost) total_cost_vre: Total annual cost of VRE in million USD (annualized cost) total_cost_capacity: Total cost of conventional + VRE capacity i.e. total_cost_new_conv + total_cost_vre total_cost_energy: Total annual cost of energy generation from all generation sources (result from economic dispatch model) total_cost: Total annual cost i.e. total_cost_capacity + total_cost_energy (annualized fixed costs of existing conventional capacity is not included here. Cancels our for additional costs of VRE) system_cost_vre: Total annual additional cost of vre i.e. (total cost of VRE scenario - total cost of zero VRE scenario)
cost_vre_pMWh: Cost of VRE per MWh of VRE after curtailment i.e. (total_cost_vre/ann_gen_vre_MWh) cost_vre_nocurt_pMWh: Cost of VRE per MWh of VRE before curtailment i.e. (total_cost_vre/ann_gen_vre_nocurt_MWh) cost_battery_pMWh: Cost of battery storage per MWh of VRE after curtailment i.e. (ann_cap_cost_battery_mil10^6/ann_gen_vre_MWh) cost_battery_nocurt_pMWh: Cost of battery storage per MWh of VRE before curtailment i.e. (ann_cap_cost_battery_mil10^6/ann_gen_vre_nocurt_MWh) cost_vre_battery_pMWh: Cost of VRE and battery storage per MWh of VRE after curtailment i.e. (cost_vre_pMWh + cost_battery_pMWh) <– USE THIS FOR EITHER OR BOTH VRE AND BATTERY cost_vre_battery_nocurt_pMWh: Cost of VRE and battery storage per MWh of VRE before curtailment i.e. (cost_vre_nocurt_pMWh + cost_battery_nocurt_pMWh) <– USE THIS FOR EITHER OR BOTH VRE AND BATTERY
capacity_value_pMWh: Capacity value of VRE target per VRE generation after curtailment i.e. difference between annualized costs of new conventional capacity of VRE scenario and zero VRE scenario per MWh of VRE generation energy_value_pMWh: Energy value of VRE target per VRE generation after curtailment i.e. difference between total annual energy costs of VRE scenario and zero VRE scenario per MWh of VRE generation total_value_pMWh: Total of capacity and energy value of VRE i.e. (capacity_value_pMWh + energy_value_pMWh)
system_cost_vre_pMWh: Additional system cost due to VRE for a VRE scenario compared to a zero VRE scenario per MWH of VRE i.e. (cost_vre_pMWh + cost_battery_pMWh - capacity_value_pMWh - energy_value_pMWh) system_cost_add_pMWh_load: Additional system cost due to VRE for a VRE scenario compared to a zero VRE scenario per MWH of Load served system_cost_add_perc: Additional overall cost of VRE as a % of no VRE system cost
total_emissions_co2_milTonnes: Total CO2 emissions in million tonnes grid_emissions_factor_kg_kWh: Grid emissions factor cost_emissions_reduction_pTonneCO2_avg: Average cost of CO2 reduction in USD per tonne of CO2 cost_emissions_reduction_pTonneCO2_marginal: Marginal cost of CO2 reduction in USD per tonne of CO2 (for every 200 GW target)
Processing the results
[1] “base”
[1] “coal_55mingen”
[1] “coal_ccgt_0mingen”
[1] “wind10LC”
[1] “wind20LC”
[1] “wind30LC”
[1] “solar10LC”
[1] “solar20LC”
[1] “solar30LC”
[1] “wind30LC_solar30LC”
[1] “wind120HH_solar1A”
[1] “battery15”
[1] “battery30”
[1] “battery15B25LC”
[1] “battery30B25LC”
[1] “battery15B50LC”
[1] “battery30B50LC”
[1] “high_cost_coal”
[1] “wind30LC_solar30LC_coalHC”
[1] “load_modD0M0_energyOnly”
[1] “load_modD50M0_energyOnly”
[1] “load_modD25M25_energyOnly”
Processed Results Table
scenario scenario_main dispatch_cost ann_gen_total_MWh
1: S0W0 base 4.607645e+12 2424553182
2: S0W200 base 3.532123e+12 2424553182
3: S50W150 base 3.583863e+12 2424553182
4: S100W100 base 3.669521e+12 2424553182
5: S150W50 base 3.817777e+12 2424553182
6: S200W0 base 4.037051e+12 2424553182
ann_gen_vre_nocurt_MWh ann_gen_vre_MWh ann_curt_vre_MWh
1: 0 0 0
2: 477166828 477093308 73516
3: 454821772 454776174 45593
4: 427217839 426206260 1011581
5: 392945562 375118489 17827089
6: 348007539 278446347 69561203
ann_gen_solarPV_MWh ann_gen_wind_MWh ann_gen_coal_MWh
1: 0 0 2226209347
2: 0 477093308 1763116977
3: 87798653 366977521 1783235771
4: 174521298 251684961 1807136233
5: 246842981 128275507 1849931247
6: 278446347 0 1945090909
ann_gen_gas_ccgt_MWh ann_gen_gas_ct_MWh ann_gen_diesel_MWh
1: 37431382 1890048 187717.2
2: 23986378 1291599 229093.1
3: 25922219 1358091 426239.2
4: 30079658 1711104 585240.6
5: 37249322 2553021 866416.0
6: 38237836 2808021 1135382.3
ann_gen_other_MWh ann_gen_hydro_MWh ann_gen_nuclear_MWh
1: 0.00 128787887 30046800
2: 1139.25 128787887 30046800
3: 0.00 128787887 30046800
4: 0.00 128787887 30046800
5: 0.00 128787887 30046800
6: 0.00 128787887 30046800
ann_discharge_bat_storage_MWh ann_charge_bat_storage_MWh
1: 0 0
2: 0 0
3: 0 0
4: 0 0
5: 0 0
6: 0 0
new_capacity_coal_MW new_capacity_gas_ccgt_MW new_capacity_gas_ct_MW
1: 151140 0 13400
2: 126720 0 28800
3: 118800 0 25000
4: 119460 0 27400
5: 126720 0 24200
6: 144540 0 11000
capacity_vre_MW capacity_solarPV_MW capacity_wind_MW
1: 0 0 0
2: 200005 0 200005
3: 200025 50015 150010
4: 200102 100079 100023
5: 200039 150010 50029
6: 200004 200004 0
capacity_bat_storage_MW energy_bat_storage_MWh scenario_build
1: 0 0 0 GW
2: 0 0 200 GW
3: 0 0 200 GW
4: 0 0 200 GW
5: 0 0 200 GW
6: 0 0 200 GW
scenario_split rank date_time scenario_econ_dispatch
1: 0-0 1 2018-01-17 00:22:52 ClcC70m
2: 0-100 2 2018-01-17 14:02:44 ClcC70m
3: 25-75 3 2018-01-17 16:51:26 ClcC70m
4: 50-50 4 2018-01-17 18:07:51 ClcC70m
5: 75-25 5 2018-01-17 19:30:53 ClcC70m
6: 100-0 6 2018-01-17 20:59:34 ClcC70m
scenario_sensitivity ann_cap_cost_coal_mil ann_cap_cost_ccgt_mil
1: ClcC70m 19006.02 0
2: ClcC70m 15935.17 0
3: ClcC70m 14939.23 0
4: ClcC70m 15022.22 0
5: ClcC70m 15935.17 0
6: ClcC70m 18176.06 0
ann_cap_cost_ct_mil ann_cap_cost_solar_mil ann_cap_cost_wind_mil
1: 873.4057 0.000 0.000
2: 1877.1705 0.000 21878.861
3: 1629.4883 3933.600 16409.829
4: 1785.9191 7871.055 10941.673
5: 1577.3446 11798.049 5472.751
6: 716.9748 15729.997 0.000
ann_cap_cost_battery_mil total_cost_new_conv total_cost_vre
1: 0 19879.42 0.00
2: 0 17812.35 21878.86
3: 0 16568.71 20343.43
4: 0 16808.14 18812.73
5: 0 17512.52 17270.80
6: 0 18893.03 15730.00
total_cost_capacity total_cost_energy total_cost system_cost_vre
1: 19879.42 70886.84 90766.26 NA
2: 39691.21 54340.36 94031.56 3265.301
3: 36912.14 55136.35 92048.49 1282.230
4: 35620.87 56454.17 92075.04 1308.773
5: 34783.32 58735.03 93518.35 2752.086
6: 34623.03 62108.48 96731.51 5965.246
cost_vre_pMWh cost_vre_nocurt_pMWh cost_battery_pMWh
1: NaN NaN NaN
2: 45.85866 45.85160 0
3: 44.73284 44.72835 0
4: 44.13996 44.03544 0
5: 46.04092 43.95214 0
6: 56.49202 45.20016 0
cost_battery_nocurt_pMWh cost_vre_battery_pMWh
1: NaN NaN
2: 0 45.85866
3: 0 44.73284
4: 0 44.13996
5: 0 46.04092
6: 0 56.49202
cost_vre_battery_nocurt_pMWh capacity_value_pMWh energy_value_pMWh
1: NaN NA NA
2: 45.85160 4.332646 34.68186
3: 44.72835 7.279860 34.63350
4: 44.03544 7.206089 33.86312
5: 43.95214 6.309745 32.39460
6: 45.20016 3.542470 31.52623
total_value_pMWh system_cost_vre_pMWh system_cost_add_pMWh_load
1: NA NA NA
2: 39.01451 6.844155 1.3467638
3: 41.91336 2.819476 0.5288523
4: 41.06921 3.070751 0.5397998
5: 38.70434 7.336577 1.1350898
6: 35.06870 21.423322 2.4603485
system_cost_add_perc vre_share_after_curt vre_curt vre_share_nocurt
1: NA NA NA NA
2: 0.03597483 0.1967758 0.0001540759 0.1968061
3: 0.01412673 0.1875711 0.0001002543 0.1875899
4: 0.01441916 0.1757875 0.0023678302 0.1762048
5: 0.03032058 0.1547165 0.0453677929 0.1620693
6: 0.06572096 0.1148444 0.1998841529 0.1435347
total_emissions_co2_milTonnes grid_emissions_factor_kg_kWh
1: 2065.178 0.8517769
2: 1633.143 0.6735853
3: 1652.669 0.6816387
4: 1676.757 0.6915735
5: 1719.900 0.7093680
6: 1808.244 0.7458049
total_emissions_reduced_milTonnes emissions_reduction
1: NA NA
2: 432.0350 0.2091998
3: 412.5091 0.1997450
4: 388.4218 0.1880815
5: 345.2781 0.1671904
6: 256.9348 0.1244129
cost_emissions_reduction_pTonneCO2_avg avoided_conv_cap_mw_perVREmw
1: NA NA
2: 7.557954 0.04509887
3: 3.108369 0.10368704
4: 3.369465 0.08835494
5: 7.970634 0.06808672
6: 23.216959 0.04499910
write results csv
Scenarios for individual plots
Pull scenario sensitivities
Plot themes, positions, and labels PLus choose results based on RE targets e.g. 200, 300, 400 GW OR 200, 400, 600 GW, based on user input
Create individual plots [Disabled with eval = FALSE]
Print and save the plots [Disabled with eval=FALSE]
Plot functions 1) Line plots with facets
Modify the battery plots in this code chunk
Modify the min gen plots in this code chunk
Modify the load modification plots in this code chunk